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Abstract. In a series of papers, Burdzy et. al. introduced the mirror coupling 
of reflecting Brownian motions in a smooth bounded domain D C K"* , and used 
it to prove certain properties of eigenvalues and eigenfunctions of the Neumann 
Laplaceian on D. 

In the present paper we show that the construction of the mirror couphng 
can be extended to the case when the two Brownian motions live in different 
domains Di,D2 C K'*. 

As an application of the construction, we derive a unifying proof of the 
two main results concerning the validity of Chavel's conjecture on the domain 
monotonicity of the Neumann heat kernel, due to I. Chavel (1121). respectively 
W. S. Kendall ([16]). 

1. Introduction 

The technique of couphng of reflecting Brownian motions is a useful tool used 
by several authors in connection to the study of the Neumann heat kernel of the 
corresponding domain (see 0, [3], [6], [11], [16], [l7], etc). 

In a series of paper, Krzysztof Burdzy et. al. ( [T], [5], [3], [5], [ID],) introduced 
the mirror coupling of reflecting Brownian motions in a smooth domain D C M'^ 
and used it in order to derive properties of eigenvalues and eigenfunctions of the 
Neumann Laplaceian on D. 

In the present paper, we show that the mirror coupling can be extended to the 
case when the two reflecting Brownian motions live in different domains Di , D2 C 

The main difhculty in the extending the construction of the mirror coupling 
comes from the fact that the stochastic differential equation(s) describing the mirror 
coupling has a singularity at the times when coupling occurs. In the case Di = D2 — 
D considered by Burdzy et. all. this problem is not a major problem (although 
the technical details are quite involved, see (2), since after the coupling time the 
processes move together. In the case Di 7^ D2 however, this is a major problem: 
after processes have coupled, it is possible for them to decouple (for example in 
the case when the processes are coupled and they hit the boundary of one of the 
domains). 

It is worth mentioning that the method used for proving the existence of the 
solution is new, and it relies on the additional hypothesis that the smaller domain 
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D2 (or more generally Di O D2) is a convex domain. This hypothesis allows us 
to construct an explicit set of solutions in a sequence of approximating polygonal 
domains for D2, which converge to the desired solution. 

As an application of the construction, we will derive a unifying proof of the 
two most important results on the challenging Chavel's conjecture on the domain 
monotonicity Neumann heat kernel (I12j. |16j). which also gives a possible new line 
of approach for this conjecture (note that by the results in [4, Chavel's conjecture 
does not hold in its full generality, but the additional hypotheses under which this 
conjecture holds are not known at the present moment). 

The structure of the paper is as follows: in Section [2] we briefly describe the 
construction of Burdzy et. al. of the mirror coupling in a smooth bounded domain 
domain D cM."^. 

In Section [3l in Theorem 13.11 we give the main result which shows that the 
mirror coupling can be extended to the case when D2 C Di are smooth bounded 
domains in R'' and D2 is convex (some extensions of the theorem are presented in 
Section [5]). 

Before proceeding with the proof of theorem, in Remark 13.41 we show that the 
proof of the theorem can be reduced to the case when Di = M.'^. Next, in Section 
13.11 we show that in the case D2 = (0,oo) C Di — R the solution is essentially 
given by Tanaka's formula (Remark 13. 5p . and then we give the proof of the main 
theorem in the 1-dimensional case (Proposition 13. 6| ). 

In Section 13. 2[ we first prove the existence of the mirror coupling in the case 
when D2 is a half-space in R'' and Di = W'' (Lemma 13. 8p . and then we use this 
result in order to prove the existence of the mirror coupling in the case when D2 
is a convex polygonal domain in M'' and Di = (Theorem 13. 9p . Some of the 
properties of coupling, essential for the extension to the general case are detailed 
in Proposition l3.10l 

In Section 2] we give the proof of the main Theorem 13. II The idea of the proof is 
to construct a sequence (F^", Xt) of mirror couplings in [Dn, Di), where Dn D2 
is a sequence of convex polygonal domains in R*^, and to use the properties of the 
mirror couplings in polygonal domains (Proposition I3.10p in order to show that 
the sequence Y^" converges to a process Yt, which gives the desired solution to the 
problem. 

The last section of the paper (Section[S|) is devoted to discussing the applications 
and the extensions of the mirror coupling constructed in Theorem 13.11 First, in 
Theorem 15.31 we use the mirror coupling in order to give a simple, unifying proof 
of the results of I. Chavel and W. S. Kendall on the domain monotonicity of the 
Neumann heat kernel (Chavel's Coniecture lS.ip . The proof is probabilistic in spirit, 
relying on the geometric properties of the mirror coupling. 

In Remark 15.51 we discuss the equivalent analytic counterpart of the proof The- 
orem [531 which might give a possible new line of approach for extending Chavel's 
conjecture to other classes of domains. 

Without giving all the technical details, we discuss the extension of the mirror 
coupling to other classes of domains (smooth bounded domains Di,2 C R'' with 
non-tangential boundaries, such that Di H D2 is a convex domain). 

The paper concludes with a discussion on the (non) uniqueness of the mirror 
coupling. It is shown here that the lack of uniqueness is due to the fact that after 



MIRROR COUPLING OF RBM AND APPLICATIONS 



3 



coupling, the processes might decouple, not only on the boundary of the domain, 
but even when they are inside of it. 

The two basic solutions give rise to the sticky, respectively non-sticky mirror 
couplings, and there is a whole range of intermediate possibilities. The stickiness 
refers to the fact that after coupling the processes "stick" to each other as long as 
possible (this is the coupling constructed in Theorem 13. 1|) hence the name "sticky" 
mirror coupling, or they can split apart immediately after coupling, in the case of 
"non-sticky" mirror coupling, the general case (weak/mild mirror couplings) being 
a mixture of these two basic behaviors. 

We developed the extension of the mirror coupling having in mind the application 
to Chavel's conjecture, for which the sticky mirror coupling is the "right" tool, but 
perhaps the other mirror couplings (the non-sticky and the mild mirror couplings) 
might prove useful in other applications. 



2. Mirror couplings of reflecting Brownian motions 

Reflecting Brownian motion in a smooth domain D C M"^ can be defined as a 
solution of the stochastic differential equation 



(2.1) Xt=x + Bt + f VDiXs)dL 

Jo 



X 

s ' 



where Bt is a d-dimensional Brownian motion, vo is the inward unit normal vector 
field on dD and is the boundary local time of Xt (the continuous non-decreasing 
process which increases only when Xt G dD). 

In (Tj , the authors introduced the mirror coupling of reflecting Brownian motion 
in a smooth domain D C K"^ (piecewise domain in with a finite number of 
convex corners or a domain in R"^, d>3). 

They considered the following system of stochastic differential equations: 



(2.2) Xt = x + Wt+ VD{Xs)dL 



t 

X 



(2.3) Yt = y + Zt+ UD{Xs)dL 



Y 



(2.4) Zt = Wt-2 [ ^\ {Y, - Xs) ■ dWs 

Jo \\Ys - Xs\\ 

for t < where ^ = inf {s > : Xg = Yg} is the coupling time of the processes, 
after which the processes X and Y evolve together, i.e. Xt = Yt ioi t > ^ and 
Zt^Z^ + lt>^ {Wt - W^). 

In the notation of [1 , considering the Skorokhod map T : C([0, oo) : R'^) — >■ 
C ([0, oo) : D), we have X = r (x + W) and Y = r {y + Z), and the above system 
reduces to 

(2.5) Zt= G{r{y + Z)^-r{x + W)J dWg + lt>^ {Wt -w^), 

Jo 

where ^ = inf {t > -.T {x + Ws) =r {y + Zg)}, for which the authors proved the 
pathwise uniqueness and the strong uniqueness of the process Zt (given the Brow- 
nian motion Wt). 
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In the above C? : M — >■ A^dxd denotes the function defined by 

(2.6) G(z) = |^(t)' if^^O 

[ 0, ifz = 

where for a unitary vector rn E M"*, H (m) represents the Unear transformation 
given by the d x d matrix 

(2.7) H (m) = I -2m m', 
that is 

(2.8) H (m) V ~ V — 2 {m ■ v) m 

is the mirror image of w G R*^ with respect to the hyperplane through the ori- 
gin perpendicular to m (m' denotes the transpose of the vector m, vectors being 
considered as column vectors). 

The pair {Xt^Yt)^^^^ constructed above is called a mirror coupling of reflecting 
Brownian motions in D starting at a;, j/ G D. 

Remark 2.1. The relation (|2.4I) can be written 

which shows that for i < ^ the increments of Zt are mirror images of the increments 
of Wt with respect to the line of symmetry Mt of Xt and Yt, which justifies the 
name of mirror coupling. 



3. Extension of the mirror coupling 

The main contribution of the author is the observation that the mirror coupling 
introduced above can be extended to the case when the two reflecting Brownian 
motion have different state spaces, that is when Xt is a reflecting Brownian motion 
in Di and Yt is a reflecting Brownian motion in D2 ■ Although the construction can 
be carried out in a more general setup (see the concluding remarks in Section [SJ, in 
the present section we will restrict to the case when one of the domains is strictly 
contained in the other one. 

The main result is the following: 

Theorem 3.1. Let -D1.2 C K'^ be smooth bounded domains (piecewise -smooth 
boundary with convex corners in M^, or -smooth boundary in , d > 3 will 
suffice) with D2 C Di and D2 convex domain, and let x G Di and y G D2 be 
arbitrarily fixed points. Given a d-dimensional Brownian motion (M^t)t>o starting 
at on a probability space {^,J-, P), there exists a strong solution of the following 
system of stochastic differential equations 

(3.1) Xt = x + Wt+ f vmiXs) dLf 

Jo 

(3.2) Yt = y + Zt+ f iyDAYs)dLj 

(3.3) Zt = f G{Y,-X,)dWs 

Jo 



MIRROR COUPLING OF RBM AND APPLICATIONS 



5 



or equivalent 

(3.4) Zt^ j^G(v{y + Z)^ - T (x + W))j dWs, 

where T and T denote the corresponding Skorokhod maps which define the reflecting 
Brownian motion X = T{x + W) in Di, respectively Y = T [y ~\- Z) in D2, and 
G ^ M.dxd denotes the following modification of the function G defined in the 
previous section: 




ifz^Q 
ifz = 



Remark 3.2. As it will follow from the proof of the theorem, with the choice of 
G above, in the case Di = D2 = D the solution of the equation p.4p given by 
the theorem above is the same as the solution of the equation (12.51) considered by 
the authors in [1] (as pointed out by the authors, the choice of G (0) is irrelevant 
in this case), and therefore the above theorem is a natural generalization of their 
result to the case when the two processes live in different spaces. We will refer 
to a solution Xt,Yt given by the above theorem as a mirror coupling of reflecting 
Brownian motions in Di, respectively D2, starting from {x,y) € Di x D2 with 
driving Brownian motion Wt- 

As we will see in Section[5l without additional assumptions, the solution of (|3.4p 
is not pathwise unique. This is due to the fact that the stochastic differential 
equation has a singularity at the origin (i.e. at times when the coupling occurs); 
the general mirror coupling can be thought as depending on a parameter which 
is a measure of the stickiness of the coupling: once the processes Xt and Yt have 
coupled, they can either move together until one of them hits the boundary {sticky 
mirror coupling - this is in fact the solution constructed in the above theorem), or 
they can immediately split after coupling (non-sticky mirror coupling) , and there is 
a whole range of intermediate possibilities (see the discussion at the end of Section 

ED. 

As an application, in Section [5] we will use the former mirror coupling (the 
sticky mirror coupling) to give a unifying proof of Chavel's conjecture on the do- 
main monotonicity of the Neumann heat kernel for domains -Di,2 satisfying the ball 
condition, although the other possible choices for the mirror coupling might prove 
useful in other contexts. 

Before carrying out the proof, we begin with some preliminary remarks which 
will allow us to reduce the proof of the above theorem to the case Di = R'^. 

Remark 3.3. The main difference from the case when Di = D2 = D considered by 
the authors in [1] is that after the coupling time ^ the processes Xt and Yt may 
decouple. For example, if t > ^ is a time when Xt = Yt € dD2, the process Yt 
being conditioned to stay in D2 , receives a push in the direction of the inward unit 
normal to the boundary of D2, while the process Xt behaves like a free Brownian 
motion near this point (we assumed that D2 is strictly contained in Di), and 
therefore the processes X and Y will drift apart, that is they will decouple. Also, 
as shown in Section O because the function G has a discontinuity at the origin, 
it is possible that the solutions decouple even inside the domain D2, so, without 
additional assumptions, the mirror coupling is not uniquely determined (there is 
no pathwise uniqueness of (|3.4p ). 
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Remark 3.4. To fix ideas, for an arbitrarily fixed e > chosen small enough such 
that £ < dist 91)2), we consider the sequence (Cn)r!>i '^^ coupling times and 

the sequence (t„)^>q of times when the processes are e-decoupled (e- decoupling 
times, or simply decoupling times by an abuse of language) defined inductively by 

^„ = inf {t > Tn-i : Xt ^ Yt} , 
Tn = mf{t> U-\Xt-Yt\> e}, 

where tq — and = C is the first coupling time. 

To construct the general mirror coupling (that is, to prove the existence of a 
solution to p.l|) - p.3p above, or equivalent to (13. 4p ). we proceed as follows. 

First note that on the time interval [0,^], the arguments used in the proof of 
Theorem 2 in [1] (pathwise uniqueness and the existence of a strong solution Z of 
(1231)) do not rely on the fact that Di = D2, hence the same arguments can be used 
to prove the existence of a strong solution of p.4p on the time interval [0, ^1] = [0, £]. 
Indeed, given Wt, p.ip has a strong solution which is pathwise unique (the reflecting 
Brownian motion Xt in -Di), and therefore the proof of pathwise uniqueness and 
the existence of a strong solution of (13.41) is the same as in 1 considering D = D2- 
Also note that as pointed by the authors, the value G (0) is irrelevant in their 
proof, since the problem is constructing the processes until they meet, that is for 
Yt — Xt^ 0, for which the definition of G coincides with p.5|) . 

Next, assuming the existence of a strong solution to p.4p on [Ci , ti] (and therefore 
on [0,Ti]), since at time ti the processes are £ > units apart, we can apply again 
the results in [1] (with Bt = Bt+n — and starting points Xt-^ and ) in order 
to obtain a strong solution of (|3.4p on the time interval [ti,^2], and therefore by 
patching we obtain the existence of a strong solution of p.4p on the time interval 

[0,6]. 

For an arbitrarily fixed t > 0, since only a finite number of coupling/decoupling 
times ^„ and r„ can occur in the time interval [0,t] (we use here the fact that D2 
is strictly contained in Di), it follows that there exists a strong solution to (|3.4p on 
[0,t] for any t > (and therefore on [0,cx))), provided we show the existence of a 
strong solution of (|3.4p on , t„] , n> 1 . 

In order to prove this claim, it suffices therefore to show that for any starting 
points X = y £ D2 of the mirror coupling, there exists a strong solution to (|3.4p 
until the £-decoupling time n. Since e < dist {dDi,dD2), it follows that the process 
Xt cannot reach the boundary dDi before the £-decoupling time ri, and therefore 
we can consider that Xt is a free Brownian motion in R'', that is we can reduce the 
proof of Theorem 13. II to the case when Di = R'^. 

We will first give the proof of the in the 1-dimcnsional case, then we will extend 
the construction to polygonal domains in R'', and we will conclude with the proof 
in the general case. 



3.1. The 1-dimensional case. From Remark 13.41 it follows that in order to con- 
struct the mirror coupling in the 1-dimensional case, it suffices to consider Di — R 
and £>2 = (0, a) and to construct a strong solution for t < n = inf {s > : |Xs — | > 
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of the following system: 

(3.6) Xt = x + Wt 

(3.7) Yt = x + Zt + L( 



(3.8) 



Zt = I G{Ys- X,) dWs 
Jo 



for an arbitrary choice x € [0, a] of the starting point of the mirror coupling, where 
e € (0,a) is arbitrarily small, (W^t)4>o ^ 1-dimensional Brownian motion starting 
at Wq = and the function G:R— >-A^ixi=IRis given in this case by 



Gix) 



-1, if X 7^0 

+1, ifx = 



Remark 3.5. Before proceeding with the proof, it is worth mentioning that the 
heart of the construction is Tanaka's formula. To see this, consider for the moment 
a = oo, and note that Tanaka formula 

\x + Wt\=x+ f sgn{x + Ws)dWs + L"tix + W) 
Jo 

gives a representation of the reflecting Brownian motion |a; + in which the 
increments of the martingale part of |a; + Wt\ are the increments of Wt when x + 
Wt G [0, oo), respectively the opposite (minus) of the increments of Wt in the 
opposite case (L^ (x + W) denotes here the local time at of x + Wt). 

Noting that the condition x + Wt & [0, oo) is the same as \x + Wt\ = x + Wt, 
from the definition of the function G it follows that the above can be written in the 
form 

\x + Wt\^x+ f Gi\x + Ws\-ix + Ws))dWs+L^+^, 
Jo 

which shows that a strong solution to (|3.6I) - p.Sp above (in the case a — oo) is 
given exphcitly hy Xt = x + Wt and Yt — \x + Wt\ (and Zt — /p sgn {x + Ws) dWs). 

We have the following: 

Proposition 3.6. Given a l-dimensional Brownian motion {Wt)i^Q starting at 
Wo — 0, a strong solution to iS. 6]) - \3. f^) above fort < n = inf {s > : \Xs — Ys\ > e} 
is given by 

Xt = x + Wt 
Yt = \a-\x + Wt-a\\ 
Zt = /o sgn {Ws) sgn (a - W,) dW, 

where 

-1, tfx>0 



sgn {x) ^ 



-1, ifx<0 



Proof. Since £ < a, it follows that for t < n we have Xt = x + Wt & {—a, 2a), and 
therefore 

-{x + Wt), x + Wte i-a, 0) 
(3.9) Yt = \a-\x + Wt-a\\^ { x + Wt, x + Wte[0,a] 

2a-x-Wt, x + Wte{a,2a) 
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Applying the Tanaka-Ito formula to the function / (z) = |a — |z — a|| and to the 
Brownian motion Xt = x + Wt, for t < ti we obtain 

Yt = x+ f sgn {x + Ws)sgn (a- x-Ws)d{x + Ws) + L"t-L'} 
Jo 

= x+ [ sgn{x + Ws)sgnia-x-Ws)dWs+ [ vdAYs) d{Ll + Li) , 



where = sups<j (x + Ws)^ and ~ sup^<( {x + Ws — a)^ are the local times 
oi X + Wt at 0, respectively at a, and vd^ (0) = +1, vd^ [o) = ^1- 
From (13.91) and the definition of G we have 



-1, x + W", e(-a,0) 
sgn{x + Ws)sgn{a-x-Ws) = { +1, x + Ws^[Q,a] 

-1, x + Ws& (a, 2a) 

/ +1, Xs=Ys 
\ -1, Xs^Ys 

= G{Y,~X,), 
and therefore the previous formula can be written equivalently 



where 



Yt = x + Zt+ I VDAYs)dL 
Jo 

Zt^ [ G{Ys~Xs)dWs 
Jo 



and LY = + LI is a. continuous nondecreasing process which increases only when 
X + Wt {0, a}, that is only when Yt G dD2, which concludes the proof. □ 

3.2. The case of polygonal domains. In this section we will consider the case 
when Z?2 C Di C M'' are convex polygonal domains (convex domains bounded by 
hyperplanes in R'^). From Remark 13.41 it follows that we can consider Di = M.'^ and 
therefore it suffices to prove the existence of a strong solution to 



(3.10) Xt = Xo + Wt 



(3.11) Yt - Yo + Zt+ f iyDAY,)dLY 

Jo 

(3.12) Zt ^ I G{Ys-Xs)dWs 

Jo 

or equivalently 

(3.13) Jo^ + " ^" " ^0 

where Wt is a d-dimensional Brownian motion starting at Wo — and Xq = Yo G 

D2. 

The construction relies on the skew product representation of Brownian motion 
in spherical coordinates, that is 

(3.14) Xt^RtQa,, 
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where Rt ^ \Xt\ £ BES (d) is a Bessel process of order d and 9f e EM (5'* ^) is 
an independent Brownian motion on the unit sphere S'''^^ in R"*, run at speed 







(3.15) at 
which depends only on Rt- 

Remark 3.7. One way to construct the Brownian motion Qt — Qt ~^ on the unit 
sphere S"^^^ C M'* is to proceed inductively on c? > 2, using the skew product 



representation of Brownian motion on the sphere Qf ^ G S"' 

ef-i = (cos singlet') 



d-1 



[see \T5\ ) 



where 9^ e LEG (d — 1) is a Legendre process of order d — 1 on [0, n], and 8f ^ £ 
S'^^'^ is an independent Brownian motion on S'^^'^, run at speed 



at 



1 



sm- 



ds. 



Therefore, considering independent processes 9} , . . . 9f~^, where 9^ e LEG [d — i) 
on [0,7r] (i = Ijd — 2) and 9f~^ a 1-dimensional Brownian (9( — {cos 9], sva.9l) S 
is a Brownian motion on S^), we have 



e 



— ^Luo , Dili yj^ i/^ , Sill Sill lud 

or equivalent, in spherical coordinates, B^^"*^ € S^~^ is given by 

(3.16) e^l = (0^...,^?^^0^l). 

To construct the solution we first consider first the case when D2 is a half-space 
%+ = {{z^^...,z'^) cW^:z'^> 0}. 

Given an angle S M, we introduce the rotation matrix R{ip) G Mdxd which 
leaves invariant the first d—2 coordinates and rotates clockwise by the angle a the 
remaining 2 coordinates, that is 

\ 



, sm ut ■ ■ ■ sm 



nd-l , 



(3.17) 



R{a) 



( 1 




V 















1 

COS Lp — sin (/3 
sin COS j 



We have the following: 
Lemma 3.8. Let D2 = = { (z^, . . . , z"*) e M'' : z'* > 0} and assume that 
(3.18) Ya = R{ipo)Xo 

for some (po G M. 

Consider the reflecting Brownian motion 9f^^ on [0, tt] with driving Brownian 
motion 9^~^ , where 9f~^ is the (d — 1) spherical coordinate ofG {Yq — Xq) Xt, given 
by [3l4\ l ~ l[3A6\) above, that is: 

9f-^ = 6lf-i + L° [9'^-^^ - LJ (9'^-^^ , t>0, 
and Lt ^f)"^^^^ , Lt ^f)"^^^^ represent the local times of 9"^^^ at 0, respectively at tt. 
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A strong solution to i3.10\) - i3.12\) above is explicitly given by 
(3.19) y^^l Ri^.)GiY,-X,)X„ t<C 

where ^ = inf {t > : Xt — Yt} is the coupling time, the rotation angle ipt is given 
by 



and for z = (z^, . . . , z'') G M'^ we denoted by |z|^ = (^z^ , z^, . . . , \z'^\) . 




Figure 1 . The mirror couphng of a free Brownian motion Xt and 
a reflecting Brownian motion Yt in the half-space "HJ". 



Proof. Recall that for to e M''— {0}, G (to) v denotes the mirror image (symmetric) 
of ?; G M.'^ in the hyperplane through the origin, perpendicular to m. 

By Ito formula, we have 
(3.20) 

Yt^^ = Yo+J^ R{^s)G{Yo~Xo)dXs + R (^^s + ^) G {Yo - Xo) dL^ 

Note that the composition Ro G {a symmetry followed by a rotation) is a sym- 
metry, and since \Yt\ = \Xt\ for all t > 0, it follows that Xt and 1* are symmetric 
with respect to a hyperplane passing through the origin for all t < ^; therefore, 
from the definition (|3.5p of G it follows that we have Yt = G {Yt — Xt) Xt for all 

Also note that when [6'^^^ ) increases, Y^ € dD2 and we have 



R (^ifis + I) G (Yo - Xo) X,=R (I) Ys = i^d, {Ys 
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Y 
s > 



and if [9 j increases, Yg G 9I?2 and we have 

R + I) G (Fo - ^o) Xs^R (I) n = -VD, (n) , 
It follows that the relation (I3.20p above can be written in the form 

Ytr.i = Yo+ G{Ys~ Xs) dXs + / i^D, {Ys) dL 
Jo Jo 

where Lf — Lj ^^''^^^ + (j)"^^^^ is a continuous non-decreasing process which 

increases only when Yt E dD2 , and therefore Yt given by p.l9p is a strong solution 
to ((XTU)) - (|XT^ for t < ^. 

For t > ^, we have Yt = \Xt\^ = {X^ , Xf , . . . ,\X^\) , and similar to the 1- 
dimensional case, by Tanaka formula we obtain: 

(3.2]y),v« = n + y (l,...,l,sgn(xf))dX, + y (0,0,...,1)L?(X'') 

= n + y G{Ys-Xs)dXs + J vdAYs)LJ, 



= (l,...,l,sgn(X,^)). 

LY — Lf {X'^^ in (|3.2ip is a continuous non-decreasing process which increases 
only when Yt € 9I?2 (we denoted by {X'''^ the local time at of the last cartesian 
coordinate X'^ of X), which shows that Yt also solves p. 101) - p.l2p for t > ^, and 
therefore Yt is a strong solution to (|3.10p - p.l2p for t > 0, concluding the proof. □ 

Consider now the case of a general polygonal domain D2 C M'^. We will show 
that a strong solution to (|3.10p - (I3.12p can be constructed from the previous lemma, 
by choosing the appropriate coordinate system. 

Consider the times {o'n)n>o which the solution Yj hits different bounding 
hyperplanes of dD2, that is do = inf {s > : G dD2} and inductively 

Yt e 9I?2 and Yt,Y(7„ belong to different^ 



(3.22) (Tn+i inf'^i (t„ . bounding hyperplanes of 9Z?2 

If Xq = Yq E dD2 belong to a certain bounding hyperplane of D2 , we can chose 
the coordinate system so that this hyperplane is Hd = {{^^ , ■ ■ ■ , z'^) G M'' : z'^ = O} 
and D2 C "Hj", and we let Hd be any bounding hyperplane of D2 otherwise. 

Then, on the time interval [ao, cti), the strong solution to (I3.10p - (|3.12p is given 
explicitly by (|3.19p in Lemma 13.81 

If fJi < 00, we distinguish two cases: Xa--^ — Y^^ and X^,-^ ^ Yo-i- Let % denote 
the bounding hyperplane of D which contains Yo-i, and let v-u denote the unit 
normal to T-L pointing inside D2. 

If Xo-i = ^0-1 G choosing again the coordinate system conveniently, we may 
assume that H is the hyperplane is Tid — {(z^, . . . , z'^) : z'^ — Q^, and on 
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the time interval [cri,cr2) the couphng {Xai+t,Ya^+t)t£[o,a2~ai) given again by 
Lemma 13.81 

If Xo-i ^ Yo-i € Ji, in order to apply the lemma, we have to show that we can 
choose the coordinate system so that the condition (|3.18p holds. If Yo-i — X^-^ is a 
vector perpendicular to "H, by choosing the coordinate system so that H — Hd = 
{(2;^, . . . , z"^) G M'^ : z"* = 0}, the problem reduces to the 1-dimensional case (the 
first d — 1 coordinates of X and Y are the same), and it can be handled as in 
Proposition 13.61 bv the Tanaka formula. The proof being similar, we omit it. 

If ^ Yfjx G and F^i —X^^ is not orthogonal to consider Xa^ = pr^ X„^ 
the projection of X^^ onto and therefore X^^ ^ Y^^. The plane of symmetry 
of Xci and Ycj^ intersects the line determined by X^^ and Y^^ at a point, and we 
consider this point as the origin of the coordinate system (note that the intersection 
cannot be empty, for otherwise the vectors Y^^ — X„^ and Y„^ — X^^ were parallel, 
which is impossible since then Y^i — X^^ , Y^^ — X„^ and Y„^ — X^,^ , X^,^ — X^^ were 
perpendicular pairs of vectors, contradicting X„^ 7^ Y^i - see Figure [2]). 




Figure 2. Construction of an appropriate coordinate system. 

Choose a orthonormal basis {ei, . . . , ed\ in I^'' such that — v-h is the normal 
vector to "H pointing inside D2, &d-\ — — — |" (^o-i ~ ^ax) is a unit vector 

lying in the 2-dimensional plane determined by the origin and the vectors and 
Yo-j — Xf,^^ and {ei, . . . , ed-2} is a completion of {e^-i, ed} to a orthonormal basis 
in R'^ (see Figured]). 

Note that by the construction, the vectors ei, . . . , ed-2 are orthogonal to the 
2-dimensional hyperplane containing the origin and the points X^-^ and Y^-^ , and 
therefore Xa-^ and Y^-^ have the same (zero) first d — 2 coordinates; also, since Xa--^ 
and Yo-i are at the same distance from the origin, it follows that Yo-i can be obtained 
from Xcr^ by a rotation which leaves invariant the first d — 2 coordinates, which 
shows that the condition (|3.18p of Lemma [3^ is satisfied. 

Since by construction the bounding hyperplane H of D2 at is given by 
-Hd €R'^:z'^ = 0} and D2 C = { (z^, . . . , z'*) €R'^:z'^> O}, 

we can apply Lemma 13.81 and deduce that on the time interval [ai,a2) a solution 
to - is given by (X,,+t, r,,+t)te[o,..~ai)- 

Repeating the above argument we can construct inductively (in appropriate co- 
ordinate systems) the solution to (|3.10p - ()3.12p on any time interval [(t„, cr„_|.i), 
n > 1, therefore obtaining a strong solution to (|3.10p - (|3.12p defined for t >0. 
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We summarize the above discussion in the following: 

Theorem 3.9. If D2 C is a convex polygonal domain, for any = Iq G D2, 
there exists a strong solution to \3. 1 0|) - 13.12\} above. 

Moreover, between successive hits of different bounding hyperplanes of D2 (i.e. 
on each time interval [(7„,(t„+i) in the notation above) and for an appropriately 
chosen coordinate system, the solution is given by Lemma\3.- 



We will refer to the solution (Xt,Yt)j^g constructed in the previous theorem as 
a mirror coupling of reflecting Brownian motions in {W^ , D2) with starting point 

If Xf Yt, the hyperplane Mt of symmetry between Xt and Yt , passing through 
with normal mt = \Yt-Xt\ ~ Xt), will be referred to as the mirror of the 
coupling. For definiteness, when Xt — Yt we let Mt denote any hyperplane passing 
through Xt = Yt, for example we choose Mt such that it is a left continuous function 
with respect to t. 

Some of the properties of the mirror coupling are contained in the following: 

Proposition 3.10. If D2 C is a convex polygonal domain, for any Xq — Yq ^ 

D2, the mirror coupling given by the previous theorem has the following properties: 

i) // the reflection takes place in the boundary hyperplane % of D2 with inward 
unitary normal V'n, then the angle /l(mt;h'-H) decreases monotonically to 
zero. 

ii) When processes are not coupled, the mirror Mt lies outside D2 . 

iii) Coupling can take place precisely when Xt G dD2. Moreover, if Xt € D2, 
then Xt ^Yt. 

iv) If Da C Dp are two polygonal domains and (Yt";Xt), {Yt^;Xt) are the 
corresponding mirror coupling starting from x G Da, for any t > we have 

(3.23) sup |r" -Yf\< Dist (D°',D^) := max \xa -xb\. 

s<t x^edDa,,XfiGdDfi ^' 



Proof, i) In the notation of Theorem 13.91 on the time interval [aQ,ai) we have 
Yt — Xt so /I {mt, v-h) = Oj and the claim is verified in this case. 

On an arbitrary time interval [cr„,(T„+i), in an appropriately chosen coordinate 
system, Yt is given by Lemma l3^ For t < (^,Yt is given by the rotation R{<pt) of 
G (Fo — -'^o) Xt which leaves invariant the first {d — 2) coordinates, and therefore 

Z {mt,vn) = ^ {mo, vn) + — 

which proves the claim in this case (note that before the coupling time ^ only one 
of the non-decreasing processes and Lt is not identically zero) . 

Since for t > ^ we have Yt = {X}, \Xt\^, we have Z {mt,vn) = which 
concludes the proof of the claim. 

ii) On the time interval [(Tq, (Ti) the processes are coupled, so there is nothing to 
prove in this case. 

On the time interval [ai,a2), in an appropriately chosen coordinate system we 
have Yt = (X^, \Xf\), thus the mirror Mt coincides with the boundary hy- 
perplane Hd — { (2^! ■ ■ • , z"^) € M'' : z'^ = 0} of D2 where the reflection takes place, 
thus MtC] D2 — in this case. 



14 



MIHAI N. PASCU 



Inductively, assume the claim is true for t < (t„ . By continuity, Mcr„ H D2 = , 
thus D2 lies on one side of Mcr„. By the previous proof, the angle Z{mt,i'n) 
between rrit and the inward unit normal to bounding hyperplane H of D2 where 
the reflection takes place decreases to zero; since D2 is a convex domain, it follows 
that on the time interval [(T„, ctk+i) we have Mt D D2 = 0, concluding the proof. 

iii) The first part of the claim follows from the previous proof (when the processes 
are not coupled, the mirror (hence Xt) lies outside D2', by continuity, it follows that 
at the coupling time ^ we must have = E dD2)- 

To prove the second part of the claim, consider an arbitrary time interval [(t„, (Jn+i) 
between two successive hits of Yt to different bounding hyperplanes of D2. In an 
appropriately chosen coordinate system, Yt is given by Lemma 13.81 After the cou- 
pling time ^, Yt is given by Yt = {X} , . . . , |-'S^t'|), and therefore if Xt € D2 (thus 
Xf > 0) we have Yt = {X^,..., Xf) = Xt, concluding the proof. 

iv) Let Af" and Aff denote the mirrors of the coupling in D", respectively , 
with the same driving Brownian motion Xt . 

Since Y" and Xt are symmetric with respect to A/", and Yf and Xt are sym- 
metric with respect to Aff , it follows that Yf is obtained from Yf by a rotation 
which leaves invariant the hyperplane Af" n Aff , or by a translation by a vector 
orthogonal to Af" (in the case when and Aff are parallel). 

The angle of rotation (respectively the vector of translation) is altered only when 
either Y^°' or Yf are on the boundary of Da, respectively Dp. Since Da C Dfj 
are convex domains, the angle of rotation (respectively the vector of translation) 

decreases when e Dp or when Y^°' e dDa and {vf - ■ vd^ {Y^°') > 
(in these cases Yf and Ft" receive a push such that the distance 



Yt'' - Yl 



IS 



decreased), thus the maximum distance 



- Yl 



is attained when Yl" € dDa 



and {yI - Fj") • vd^ (F") < 0, and the formula follows. □ 

4. The proof of Theorem 13.11 

By Remark |3.4[ it suffices to consider the case when Di ~ and D2 C M'^ is a 
convex bounded domain with smooth boundary. To simplify the notation, we will 
drop the index and write D for D2 in the sequel. 

Let {Dk)j.^i be an increasing sequence of convex polygonal domains in R'' with 

D71 C Dn+i and U„>i Ai ^ D. 

Consider {Yl^,Xt)^yQ a sequence of mirror couplings in (£)„,R'^) with starting 
point X G Di, with driving Brownian motion (W^t)j>Q, Wq = given by Theorem 

1311 

By Proposition 13. 10[ for any t > we have 

sup |y™ - < Dist {Dn, An) - max |s„ - 0, 

s<t x„€dD„.Xm&dDm n.m— >oo 

hence Yj" converges a.s. in the uniform topology to a continuous process Yt. 

Since (Y"")„>-^ are reflecting Brownian motions in (£'«)„>! and £>„ D, the law 
of Yt is that of a reflecting Brownian motion in D, that Yt is a reflecting Brownian 
motion in D starting a,t x E D (see |BJ). Also note that since Y/^ are adapted to 
the filtration — {J^t)t>o generated by the Brownian motion Wt, so is 1*. 
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By construction, the driving Brownian motion Z" of y^" satisfies 



Jo 



t > 0. 



Consider the process 



Zt= f G{Yt-Xt)dBu 
Jo 



and note that since Y is J^^-adapted and ||G|| = 1, by Levy's characterization of 
Brownian motion, Zt is a free d-dimensional Brownian motion starting at Zg = 0, 
also adapted to the filtration . 

We will show that Z is the driving process of the reflecting Brownian motion Yj, 
i.e. we have 



Yt = x 



f G (Ys ~ Bs) dWs + Ll , t>0. 
Jo 



Note that the mapping z i — > G (z) is continuous with respect to the norm 
1 1^1 1 = ll(°^'y )ll — J2i j=i '^ij of d X d matrices at all points z £ R'^ — {0}, hence 
G (Y7 ~Xs) G (Ys - Xs) if Ys- Xs ^ 0. If Y, - X^ = 0, then either Y, = 

Bse D or yJ= Xs e dD. 

li Ys ^ Bs ^ D , since Dn -D, there exists > 1 such that Bs G Dn, hence 
Bs e Dn for aU n > iV. By Proposition [330l it follows that F," = Bs for all n>N, 
hence in this case we also have G (F/' - Bs) = G (0) G {0)^G{Ys-Bs). 

If Ys = Bs <G dD, since Dn C D we have YJ^ — Bs ^ 0, and therefore by the 
definition of G we have: 

/ ||G (y," - Xs) -G{Ys-Xs)\\^ lY,=B,edDds 

Jo 



H 



1-2- 



n" - Xs 
Yr - Xs 



1 



ds 



\Yr - Xs 



n" - Xs 

IK" - Xs 



lY,=Xs<£dDds 



r, Y^~Xs 

wr-XsW 

lY,=X^<£dDds 



y " - Xs 

\Yr ~ Xs 



lY^=x^edDds 



< 



4 / l9DiYs)ds 
Jo 



0, 



since Yt is a reflecting Brownian motion in D, and therefore it spends zero Lebesgue 
time on the boundary of D. 

Since ||G|| = 1, using the above and the bounded convergence theorem we obtain 

lim / ||G(y,"-X,)-G(n-X,)||'ds = 0, 
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and therefore by Doob's inequality it follows that 

Esup\Z^-Z/ <cE\Z^~Zt\^ <cE f \\G{YJ'^Xs)-G{Ys-Xs)\\''ds ^ 0, 

s<t Jo 

for any t > 0, which shows that converges uniformly on compact sets to Zt — 
Jl;GiYs~X,)dWs. 

From the construction, Z" is the driving Brownian motion for Y^" , that is 

Yr = x + Z^+ fvD^ {Y^dL^-, 
Jo 

and passing to the limit with n — > oo we obtain 

Yt^x + Zt + At^x+ f G{Ys-X,) dWs +At, t> 0, 

where At = lini„^oo /q ^d^ {Y^) dLj" . 

It remains to show that At is a process of bounded variation. For an arbitrary 
partition ~ to < ti < . . . ti ~ t oi [0, t] we have 



I I 
Ey\At^-At^_A - lim 


r unjY:^)dLY- 




i=l i=l 






< hmsup-E lJ" 


= limsup J ^ 


' Pd„ {s, X, y) o-„ {dy) ds 




dD„ 







where cr„ is the surface measure on dDn and the last inequality above follows from 
the estimates in [5] on the Neumann heat kernels po^ (ti v) (see the remarks 
preceding Theorem 2.1 and the proof of Theorem 2.4 in 0). 

From the above it follows that At — Yt — x — Zt \s a, continuous, J^^-adapted 
process (Yt, Zt are continuous, J^^-adapted processes) of bounded variation. 

By the uniqueness in the Doob-Meyer semimartingale decomposition of Yt - 
reflecting Brownian motion in D, it follows that 

At= I i^D{Ys)dL^, t>Q, 
Jo 

where is the local time of Y on the boundary dD, and therefore the reflecting 
Brownian motion Yt in D constructed above is a strong solution to 

Yt^x+ f G{Ys-X,)dWs+ f iyD{Ys)dL^, t>0, 
Jo Jo 

or equivalent, the driving Brownian motion Zt = G {Ys — Xs) dWs of Yt is a 
strong solution to 

Zt = J^G(T{y + Z)^ - X,) dWs, t> 0, 
concluding the proof of Theorem 13.11 
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5. Extensions and applications 

As an application of the construction of mirror coupling, we will present a uni- 
fying proof of the two most important results on Chavel's conjecture. 

It is not difficult to prove that the Dirichlet heat kernel is an increasing function 
with respect to the domain. Since for the Neumann heat kernel pu {t,x,y) of a 
smooth bounded domain 13 C M'' we have 

lim pD{t,x,y)^ \ , 

t^ca vol {D) 

the monotonicity in the case of the Neumann heat kernel should be reversed. 
The above observation was conjectured by Isaac Chavel (021), as follows: 

Conjecture 5.1 (Chavel's conjecture, [H])- Let D12 C M"^ be smooth bounded 
convex domains in M'', d > 1, and let po^ {t, x, y), PD2 {t, x, y) denote the Neumann 
heat kernels in Di , respectively D2 ■ If D2 C Di , then 

(5.1) pD^it.x.y) <pD2{t,x,y) , 
for any t >Q and x,y Cz Di . 

Remark 5.2. The smoothness assumption in the above conjecture is meant to insure 
the a.e. existence the inward unit normal to the boundaries of Di and D2, that is 
the boundary should have locally a differentiable parametrization. Requiring that 
the boundary of the domain is of class C^'" (0 < a < 1) is a convenient hypothesis 
on the smoothness of the domains Di 2- 

In order to simplify the proof, we will assume that D12 are smooth domains 
(the proof can be extended to a more general setup, by approximating D12 by less 
smooth domains). 

Among the positive results on Chavel conjecture, the most general known results 
(and perhaps the easiest to use in practice) are due to I. Chavel and W. Kendall 
(see [H], [IS]), and they show that if there exists a ball B centered at either x or 
y such that D2 C B C Di, then the inequality (|5.1I) in Chavel's conjecture holds 
true for any t > 0. 

While there are also other positive results which suggest that Chavel's conjecture 
is true (see for example [11], [M]), in [4] R. Bass and K. Burdzy showed that Chavel's 
conjecture does not hold in its full generality (i.e. without additional hypotheses). 

We note that both the proof of Chavel (the case when Di is a ball centered at 
either x or y) and Kendall (the case when D2 is a ball centered at either x or y) relies 
in an essential way that one of the domains is a ball: the first uses an integration 
by parts technique, while the later uses a coupling argument of the radial parts of 
Brownian motion, and none of them can be applied to the other case. 

Using the mirror coupling, we can derive a simple, unifying proof of these two 
important results, as follows: 

Theorem 5.3. Let D2 C Di C K'' be smooth bounded domains and assume that 
D2 is convex. If for x,y € D2 there exists a ball B centered at either x or y such 
that D2 d B d Di, then for all t > we have 

(5.2) pDAt,x,y) <pD2{t,x,y) . 

Proof. Consider x,y Cz D2 fixed and assume without loss of generality that D2 C 
B^B{y,R) C Di for some i? > 0. 
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Consider a mirror coupling {Xt,Yt) of reflecting Brownian motions in {Di,D2) 
starting at y £ D2. 

The idea of the proof is to show that at aU times Yt is at a distance from y 
smafler than (or equal) to that of Xt from y. 

To prove the claim, consider a time to > when the processes are at the same 
distance from y, that is \Ytg — y\ = \Xtg — x\. If Xtg = Ytg, for t > to the distances 
from Xt and Yt to y will remain equal until the time ti when the processes hit the 
boundary of D2, and Yt receives a push in the direction of the inward unit normal 
to the boundary of D2 ■ Since D2 is convex, this decreases the distance of Yt from 
y, and the claim follows in this case. 

If the processes are decoupled and \Ytg — y\ = \Xtg — x\, the hyperplane Mtg of 
symmetry between Xtg and Ytg passes through y, and the ball condition shows that 
we cannot have Xtg £ dDi. Therefore for t > to, the processes Xt and Yt will 
remain at the same distance from y until Yt € dD2, when the distance of Yt from 
y is again decreased by the local push received as in the previous case, concluding 
the proof of the claim. 

Therefore, for any e > we have 

P^{\Xt-y\<e)<P^{\Yt~y\<e), 

and dividing by the volume of the ball B (y, e) and passing to the limit with e \, 0, 
from the continuity of the transition density of the reflecting Brownian motion in 
the space variable we obtain 

PDi {t,x,y) < PD2 {t,x,y) , 

for any t > 0, concluding the proof of the theorem. □ 

Remark 5.4. As also pointed out by Kendall in [16j (the case when D2 is a ball), we 
note that the convexity of the larger domain Di is not needed in the above proof 
in order to derive the validity of condition (|5.ip in Chavel's conjecture. 

Remark 5.5. We also note that the above proof uses only geometric considerations 
on the relative position of the reflecting Brownian motions coupled by mirror cou- 
pling. Analytically, the above proof reduces to showing that Rt = \Xat |^ — \Ya^ |^ > 
for all < t < ^ = inf {s > : Xg = Yg}, where Rt is the solution of the following 
stochastic differential equation 

(5.3) Rt = Ro + 2 f RgdBs + 2St, 

JQ 

where Bt — J^* ^x'-y]'^ ' '^^^ 1-dimensional Brownian motion, at — A^^ is the 
inverse of the non-decreasing process At defined by 

f* 1 

At = T^ds, 

Jo \Xs-Ys\^ 

and 

St^ f ^a. • i^D^ {X^JdL^^ - [ • lyo, (Y^JdLl^. 
Jo Jo 

Perhaps a better understanding of the mirror coupling, based on the analysis of 

the local times and spent by Xt and Yt on the boundaries of Di, respectively 

D2, in connection to the geometry of the boundaries dDi and dD2 could give a 
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proof of Chavel's conjecture for some new classes of convex domains, but so far we 
were unable to implement it. 

We have chosen to carry out the construction of the mirror coupling in the case 
of smooth domains with D2 C Di and D2 convex, having in mind the application 
to Chavel's conjecture. However, although the technical details can be considerably 
longer, it is possible to construct the mirror coupling in a more general setup. 

For example, in the case when Di and D2 are disjoint domains, none of the diffi- 
culties encountered in the construction of the mirror coupling occur (the possibility 
of coupling/decoupling), so the constructions extends immediately to this case. 

The two key ingredients in our construction of the mirror coupling were the 
hypothesis D2 C Di (needed in order to reduce by a localization argument the 
construction to the case Di = M.'^) and the hypothesis on the convexity of the inner 
domain D2 (which allowed us to construct a solution of the equation of the mirror 
coupling in the case Di = R''). 

Replacing the first hypothesis by the condition that the boundaries dDi and 
9I?2 are not tangential (needed for the localization of the construction of the mirror 
coupling) and the second one by condition that Di n D2 is a convex domain, the 
arguments in the present construction can be modified in order to give rise to a 
mirror coupling of reflecting Brownian motion in {Di, 02). 



Figure 3. Generic smooth domains D1.2 C M'' for the mirror 
coupling: Di,D2 have non-tangential boundaries and Di n D2 is 
a convex domain. 

We conclude with some remarks on the non-uniqueness of the mirror coupling in 
general domains. To simplify the ideas, we will restrict to the 1-dimensional case 
when D2 = (0, 00) C Z^i M. 

Fixing X €E (0,oo) as starting point of the mirror coupling {Xt,Yt) in (Di,Z?2), 
the equations of the mirror coupling are 




(5.4) 
(5.5) 



Xt 
Yt 



x + Wt 
x + Zt + L] 



Y 



(5.6) 
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where in this case 

r -1, ifz^o 

\ +1, if z = ■ 

Until the hitting time t — {s > : Ys E dD2} of the boundary of dD2 we have 
LY = 0, and with the substitution Ut — {Yt ~ Xt), the stochastic differential 
for Yt becomes 

(5.7) Ut = f ^%^^dT4^, = f a (C/,) dWs, 

Jo ^ Jo 

where 

l-G(z) r 1, if z^O 

"(") = ^^ = | 0, if.^0 ■ 

By a result of Engelbert and Schmidt (113,) the solution of the above problem 
is not even weakly unique, for in this case the set of zeroes of the function a is 
N = {0} and cr~^ is locally integrable on M. 

In fact, more can be said about the solutions of (|5.7p in this case. It is immediate 
that both Ut = and Ut = Wt are solutions to 15.71 and it can be shown that an 
arbitrary solution can be obtained from Wt by delaying it when it reaches the origin 
(sticky Brownian motion with sticky point the origin). 

Therefore, until the hitting time r of the boundary, we obtain as solutions 

(5.8) Yt=Xt^x + Wt 
and 

(5.9) Yt = Xt-2Wt^x-Wt, 

and an intermediate range of solutions, which agree with (jS.Sp for some time, then 
switch to ((5^ (see [18;). 

Correspondingly, this gives rise to mirror couplings of reflecting Brownian mo- 
tions for which the solutions stick to each other after they have coupled (as in 
()5.8p ). or they immediately split apart after coupling (as in ()5.9|) ). and there is a 
whole range of intermediate possibilities. The first case can be referred to as sticky 
mirror coupling, the second as non-sticky mirror coupling, and the intermediate 
possibilities as weak/mild sticky mirror coupling. 

The same situation occurs in the general setup in M.'^, and it is the cause of lack 
uniqueness of the stochastic differential equations which define the mirror coupling. 
In the present paper we detailed the construction of the sticky mirror coupling, 
which we considered to be the most interesting, both from the point of view of 
constructions and of the applications, although the other types of mirror coupling 
might prove useful in other applications. 
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